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Abstract 

We investigate steady symmetric gravity water waves on finite depth. For non- 
positive vorticity it is shown that the particles display a mean forward drift, and 
for a class of waves we prove that the size of this drift is strictly increasing from 
bottom to surface. We also provide detailed information concerning the streamlines 
and the particle trajectories. This includes the case of particles within irrotational 
Jh waves. 
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This paper is concerned with the streamlines and the particle trajectories of 
j> steady gravity water waves on finite depth. Such water waves are one of the 

most common wave formations at sea. As a result of dispersion wind-generated 
gravity waves eventually sort themselves out [19,20]. Larger waves move faster 
than smaller ones and swell is generated: approximately two-dimensional wave- 
trains of periodic and symmetric waves moving with constant speed across 
the sea. The exact mathematical theory for such waves is well established, in 
particular for irrotational currents [15,26]. Those model very well the situation 
when the waves propagate into a region of still water. There are however 
experimental evidence that for some situations such a model is inadequate [24] . 
One example is tidal flow, which is more correctly modelled by waves entering 
a rotational current of constant vorticity [25]. Therefore the importance of 
water motion with a non-vanishing curl — i.e. in the presence of vorticity - 
has recently come to draw a lot of attention (see e.g. [6,10,18,30]). For us it 
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is of relevance that for arbitrary vorticity distributions there exist symmetric 
waves [7], and that any wave for which the surface profile is monotone from 
crest to trough necessarily is symmetric [3]. 

In a number of recent articles the exact behaviour of the fluid particles within 
such waves have been investigated [1,2,9,16,17]. The background is the follow- 
ing. For over a century it has been known that the very first approximation of 
steady irrotational gravity water waves display closed elliptic particle trajecto- 
ries [23] . However, as was first noted in [9] , a thorough study of the linearized 
system shows that the particle paths indeed have an elliptic shape, but are 
not closed. This is so for other types of waves too, as has been shown in [2,17]. 
Since the linearized problem can be solved explicitly, details of the particle 
paths can be more easily studied. In particular it can be seen that all particles 
traverse elliptic orbits. When vorticity is present things are not as transparent, 
not even for linear waves on a current of constant vorticity. In [13] it is e.g. 
shown that when the size of the vorticity is large, the particle paths of linear 
waves need not all be circular; some particles may move constantly forward 
along with the wave. Near the flat bed the particles however always behave 
like the classical first approximation: they move slightly forward in circular 
shapes. As will be discussed below this is all in relation to some reference 
speed, i.e. the generalized Stokes requirement (2.5). 

For exact water waves the details are far more elusive since closed expressions 
are not available. The investigations [1,16] show that even for exact irrota- 
tional water waves the particles display a mean forward drift. They also assert 
that for irrotational Stokes waves on finite as well as infinite depth all particles 
above the bed move in circular orbits. In this paper we show that the mean 
forward drift is preserved for all negative vorticity distributions. For irrota- 
tional waves and small enough rotational waves, we are able to show that this 
forward drift is strictly increasing from bottom to trough. A proof of this for 
linear waves without vorticity was given in [2]. In addition we establish some 
very nice properties of the velocity field and the particle paths, in particular so 
for irrotiational waves. We have not been able to confirm the circular orbital 
shapes for rotational waves. 

The novelty of our approach lies in the fact that via maximum principles we 
establish precise pointwise information about the velocity field and its deriva- 
tives within the entire fluid domain. In this way we extend the propositions 
in [1,16], providing further understanding even for the irrotational case. For 
steady rotational waves, to our knowledge this is the first investigation of its 
kind apart from [13]. The paper is organized as follows. Section 2 gives the 
mathematical background, while the main results are proved in Sections 3 
and 4. A synthesis and an analysis of the particle paths are given in Sec- 
tion 5, presented as two, hopefully illustrative, examples. At the end Section 6 
contains some short comments on the results. 
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2 Mathematical formulation 



Let d > be the depth below the mean water level y = 0, so that the flat 
bottom can be described by y = —d. The free surface can then be represented 
by a function 

rj e C 3 (R,R). 

We require that 7/(0) = max x( z]^{ri(x)} be the vertical coordinate of the crest 
(unique within a period). Naturally mm x& ^{r)(x)} > —d so that the trough 
is above the flat bed y = —d. In this paper we shall be concerned with the 
nontrivial case when maxr/ > minr/. The wave is steady of period L > - 
without loss of generality we may take L = 2n — and we require the surface 
profile to be monotone between crests and troughs. It is therefore symmetric 
around the crest [3], and we have that 

r)(x + 2ir) = r](x), i](x) = i](—x), and rf{x) < for x G (0, 7r). 

We let denote the fluid domain and define it as the interior of its boundary 

dfl v = {y = -d} U {x,rj(x))} xeR . 

A solution to the water wave problem is then definied as a function ip e C 2 {VL rj ) 
such that 

A^ = - 7 (V>), (x,y)en v 
\ViP\ 2 + 2gy = C, y = V (x) 

ijj = 0, y = 7]{x) 

> = -Po, y = -d, 

that is even and 27r-periodic in the x- variable. In (2.1) p is called the relative 
mass flux, the vorticity function 7 : [0, — p ] — > R is continuously differentiable, 
g > is the gravitational constant, and C is a constant related to the energy. 
The setting is that of gravitational water waves, meaning that the influence of 
capillarity is neglected in (2.1), and the water is assumed to be inviscid. The 
stream function ip is defined (up to a constant) by 

i> x = -v, i> y = u - c < 0, (2.2) 

where u, v are the horizontal and the vertical velocity, respectively, and c > 
is the constant horizontal speed of propagation. The notion of relative mass 
flux introduced in [7] captures the physical fact that the amount of water 
passing any vertical line is constant troughout the fluid domain: 

r 77(20 



/ (u(x, y) - c) dy = p , 

J — d 
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holds since u — c = ip y , and ip is constant on the surface y = r)(x) as well as 
on the bottom y = —d. 

Provided that u — c — ip y < 0, the system (2.1) can be deduced from the Euler 
equations (see e.g. [7,29] for a more detailed discussion). This assumption is 
supported by physical measurements [20]: for a wave not near breaking or 
spilling, the speed of an individual fluid particle is far less than that of the 
wave itself. For irrotational waves it is however known that there exist so called 
highest waves for which the crest is a stagnation point, i.e. Vip = 0. While the 
exact problem is still open for waves with vorticity [28], there are indications 
that for some classes of vorticity there do exist steady waves with particle 
layers not satisfying ip y < [13]. In this paper we shall however consider only 
waves that are not near breaking or stagnation, so that ip y < in Q v . 

A hodograph transform converts the free boundary problem (9) into a problem 
with a fixed boundary. Let us express the height 

h = y + d 

above the flat bed in terms of the new space variables 

q = x, p = —ip. (2.3) 
Notice that ip y < so that (2.3) is a local change of variables, with 

h = -f^. _ h = _J_ 1 

fig — — , lip — — 

iPy U — C Ipy C — U 

The above local coordinate transform is actually a global change of variables 
(see [7]) so that we can transform the problem (2.1) to these variables to obtain 

!(1 + h 2 q ) hpp - 2h p h q h pq + h 2 p h qq + -y(-p) h 3 p = in p < p < 0, 
l + h 2 q + (2gh-Q)hl = on p = 0, (2.4) 
h — on p = Po, 

with h even and of period 2n in the q variable. This is an elliptic equation (since 
hp > 0) with a nonlinear boundary condition. Instead of studying (2.1) in the 
domain Q v , which depends on 77, we investigate (2.4) in the fixed rectangle 
R = (— 7r,7r) x (po,0), looking for functions h G C 2 (R) that are 27r-periodic 
in q. Notice that knowing h(q, 0) is equivalent to knowing the free surface 
y = rj(x) as %, 0) = 77(g) + d. 

Let (x,a(x)) denote the parametrization of a (general) streamline 

{(x,y): ip(x,y) = -p}. 
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Note that since ip y < the above parametrization is sensible, and we have 



To normalize the reference frame Stokes made a now commonly accepted pro- 
posal. In the case of irrotational flow he required that the horizontal velocity 
should have a vanishing mean over a period. Stokes' definition of the wave 
speed unfortunately cannot be directly translated to waves with vorticity, a 
consequence of that dWVip ^ 0, so that ip y has different means at different 
depths in view of the Divergence theorem. In the setting of periodic waves 
with vorticity we propose the requirement 



a "Stokes' condition" at the bottom. This is consistent with deep-water waves 
(cf. [11]), and is also the choice made in [25]. Calculations performed on linear 
water waves with constant vorticity indicate that this is the natural choice, 
since that and only that choice recovers the well established bound y/gh for 
the wave speed [13]. We emphazise that (2.5) is only a convention for fixing the 
reference frame; except from the assertion of forward drift it does not change 
anything in this paper. Without such a reference it is however meaningless to 
discuss whether physical particle paths are closed or not. 



3 Streamlines and the horizontal velocity 

In this section three propositions are given. The proofs are collected at the 
end of the section. Lemma 3.1 mainly contains content known within the 
field, which we make use of in the rest of the paper. Lemma 3.2 gives detailed 
information about the horizontal velocity throughout the fluid domain. Its 
Corollary 3.3 gives a lower bound for the surface of a large class of gravity 
waves. 

Lemma 3.1 

i) Every streamline satisfies a' < for x e (0, ir), and the maximal steepness 
of the streamlines is a strictly monotone function of depth. 

ii) For 7' > and 7 < 0, the maximal horizontal velocity, maxu(x,cr(x)), 



is strictly increasing from bottom to surface. 
Hi) The vertical velocity is strictly positive for x G (0,7r), and if '7' < 0, then 





h q (Q>P)- 




(2.5) 
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The fact that h q = a' has a strict sign was first noticed in [7]. The sec- 
ond proposition is a slightly different version of the simple fact that Aip y = 
—7' (■?/>) > 0. The proof here presented is new. The requirement on 7' in 
iii) dates back to [4]. Our proof techniques rely heavily on sharp maximum 
principles. We refer the reader to the excellent sources [14,21]. 

The following lemma gives detailed information about the horizontal velocity. 

Lemma 3.2 (The horizontal velocity) 

i) If 7 < 0, then the horizontal velocity u is non-increasing from crest to 
trough, i.e. 

D x u{x,rj(x)) < for x E (0,n). (3.1) 

ii) If j = and \rf \ < l/y/3, then along any streamline (x,a(x)) holds 

D x u(x,a(x)) < for x e (0, 7r), 

and the pointwise steepness of the streamlines is everywhere decreasing 
with depth. 

iii) 7/7(0) > and 7', 7" < then 

d x u(x,y)<0 in (0, 7r) 

for waves in a neighbourhood of the bifurcation point in [8]. 

iv) If the horizontal velocity attains its maximum at the surface, then it does 
so either at the crest, or at the concave part of the surface where 

(c — u)j < g = —i]"( c ~ u ) 2 ■ 

Part i) was proved in [27], building on earlier ideas for the irrotational [26] 
and the rotational [8] cases. 

Concerning part ii) it is rather wonderful that the bound \h q \ < l/y/3 appears 
in the proof of this assertion. It is the natural number since it describes the 
interior angle at the crest of Stokes highest wave, see [28]. Once we proved 
Lemma 3.1 ii), the equality (3.5) implies that Det J7(V/i) 7^ 0, so that the 
mapping (q,p) 1— > (h q ,h p ) is everywhere locally a C 1 -diffeomorphism. The 
rotational property however seems to destroy that feature. 

Part iii) equivalently states that d y v < within the half-period (0, it), so that 
the vertical velocity behaves similarly to its linear approximation for constant 
vorticity, v(x,y) ~ sin(rr) sinh(y) (cf. [13]). 

Part iv) limits the possible classes of waves with the maximal horizontal ve- 
locity bounded away from the crest. Note also that the equality of this claim 
is based solely on the surface conditions of 2.1, and hence is unrelated to e.g. 
periodicity, symmetry, and depth. 
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Corollary 3.3 (The surface) 

If (3.1) holds, then for x G (0, 7r) we have the uniform bound rj" > ~ c „ 2 gt ? (o) > 
together with 

2 

^^-c-woy and v(x) ~ m ~ 2(c -L(o))- 

It should be pointed out that Corollary 3.3 is based solely on the surface 
conditions, hence valid for all types of gravity waves. 



Proof of Lemma 3.1 

i) For details see [7, Eq. (5.18)]. The key idea is that h q is annihilated by 



(l + hfj d 2 - 2h p h q d p d g + h 2 p d 2 + 2h q h pp d q + [&y(-p)h 2 p - 2h q h pq] 

(3.2) 

The second statement follows from applying the strong maximum prin- 
ciple to subdomains (0,7r) x (—d,p) of this half-period, 
ii) Consider h p = l/(c — u) > 0. Since h p belongs to the kernel of the 
uniformly elliptic operator 

(l + hfj d pp - 2h q h p d qp + h 2 p d qq - 2h p h qp d q + h p (2h qq + 3-fh p ) d p - 7'^, 

(3.3) 

the strong maximum principle implies that max u is never attained in the 
interior of any C 2 -subdomain of the fluid. Since 



u 



y = i ) yy = ~lf(-Po) > 



at the flat bed, it is a consequence of the Hopf boundary point lemma 
that u does not attain its maximum on the bottom. The proposition then 
follows by periodicity, 
iii) Since a 1 = —ip x /ip y and ip y < 0, the positivity of v is a direct consequence 
i). Note also that v vanishes on the flat bed. Since v G Ker{A + 7' (ip)} 
we may apply the strong maximum principle to any subdomain (0, 7r) x 
(— d, <t(x)). Symmetry yields the assertion. 



Proof of Lemma 3.2 

i) See [27][Thm 2.2] 

ii) We need to show that h qp < everywhere in q G (0,7r). Recall that 
h p = l/(c — u) and that each value of p corresponds to a streamline. 
Therefore Lemma 3.1 i) shows that h qp < along the surface for q G [0, ir]. 
At the bottom we have h qp < for q G (0, ir). This follows from Lemma i) 
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and the Hopf boundary point lemma. Along the vertical sides where q = 
and q = 7i, symmetry implies that h qp = 0. We shall establish that the 
strong maximum principle holds for h qp in q e [0, 7r], according to which 
h qp > at an interior point of the half-period would force h qp = const 
everywhere. 

We begin by differentiating (3.2) with respect to p. In that expression 
appears expressions involving h qqq and h ppp . Those are substituted using 
(3.2) and (3.3), so that on the level of third and fourth order derivatives 
of h, only such that are derivatives of h qp appears. Making repeated use of 
the elliptic equality in (2.4) it is possible to write the obtained expression 
as the following identity: 

/ 2h h 2 h \ 

(l + hfj d 2 p - 2h p h q d p d q + h 2 p d 2 q + [4h q h pp - 2h p h qp - x q + p J j d q 



, Ah 2 q h qq h p 2h pp (l + h 2 q ) \ 

+ i 3 A- 2 ^V + ^ j- )d P 

,'h 2 



| 2{h qq h pp - h 2 qp ){l-?>hl)-?>ih 



h 



2h h 2 

= yH 2 b^ h ^ + 2h ^ ~ • (3-4) 

By multiplying the elliptic equation in (2.4) with h pp and then completing 
the squares, we obtain for 7 = that 

= (h q h pp - hphqp) 2 + h 2 pp + h 2 p (h qq h pp - h 2 qp ), (3.5) 

which forces h 2 > h qq h pp everywhere in the fluid. Finally, Lemma 3.1 i) 
and the assumption guarantees that h 2 < 1/3 everywhere, 
iii) We shall determine ip xy in a manner similar to that of [12]. Since ip = 
along the surface, we have ip x = —r}ip y , and insertion into the Bernoulli 
surface condition of (2.1) yields that 

% = ~T+ff- (3 - 6) 

It follows that 

i>xy + V^yy = -d x y/(C -2gri)/(l +r] 2 ). 
On the other hand, we may differentiate ip x = r]ip y once more, obtaining 

i>xx + ^^yy + ^V^xy + i)i>y = 0- 

We now know the first derivatives of ifj in terms of n and in order to 
determine the three second derivatives we only need one more equality. 
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This is supplied by 

Vpxx + lpyy = -7(0). 

The calculation, which is tedious but easily carried out by hand, yields 
that at the surface ip xy can be determined as 

77 (2ry(C - 2g V ) + (1 - i] A ) 9 + 7 (0)y^2^ (1 + i] 2 f 2 ) 

Since < C — 2grj < C + 2gd we see that for small enough 7/, rj" the 
sign is determined by g (remember that 7(0) > by assumption), so that 
ipxy < at the surface. But ip xy obeys the maximum principle, 

(A + 7') i) xy = --y"M v > 0, 

so that a nonnegative maximum in < x < n is attained either on the 
flat bed or at the surface (since ip xy = for x = kn, k e Z). Since 
Lemma 3.1 iii) and the Hopf boundary point lemma forces that ip xy < 
for y = —d, the proof is complete, 
iv) In the spirit of [12], and according to (3.6), differentiation along the sur- 
face gives 

d, <m = mm=m^±n . (3 . 8) 

(1 + 77) 

Thus a maximum implies that either i] = or (2gi] — C)fi — g (1 + if). 
With a little manipulation, half of the result is obtained by substituting 
the second expression into the Bernoulli surface condition of (2.1). To 
obtain the inequality, first note that if the maximum is attained at the 
surface, then tp yy > at that point. But (cf. [12]) 

r,{C - 2g V ) (r) 2 - 1) + 2gif (1 + if) - 7 (0)yC^2^ (1 + iff 2 

(3.9) 

so we need only substitute (2gr) — C)rj — g (1 + if) into that expression 
to get 

^-ggT (i + VWO-Tgn - a (3 ' 10) 



Proof of Corollary 3.3 Recall (3.8). It follows from the assumption that 
^<^, waning ±Jo g (±±£)<0. (3.11) 



l + ?/ J _ C-2gr] 1 ° rfx \C-2^ 
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This can be integrated to 



i + V 2 / i 



(3.12) 



C-2g V - C-2grj(0y 



Since rj' < in (0, it) we may rearrange to obtain 



^2g( V (0)- V (x)) yJC-2gr,(Q)' 



(3.13) 



Now the assertion concerning rj is established by integrating (3.13). The bound 
on rj 1 then follows from employing the lower bound on rj to (3.12). Finally, the 
uniform bound on rj" is immediate from combining the left-hand side of (3.11) 
with (3.12). 



4 The forward drift 

In this section we prove that for waves of non-positive vorticity the forward 
drift of the fluid particles is everywhere positive within the fluid domain. For 
irrotational waves and small enough rotational waves, it is shown that the 
drift is monotone from bottom to surface. The proofs are collected at the end 
of the section. 

Theorem 4.1 

i) For 7 < there are no closed particle trajectories above the flat bed. In 
particular all fluid particles display a mean forward drift. 

ii) 7/7 = with \r)'\ < l/y/3 then the mean forward drift is strictly increasing 
from bed to surface. 

Hi) If 7 < then for all waves in a neighbourhood of the bifurcation point 
found in [7], the mean forward drift is strictly increasing from bed to 
surface. 

Lemma 4.2 For (x,a(x)) a non-trivial streamline, and 7 < 0, the quantity 




is non-decreasing with depth, and 




(4.1) 
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Proof of Lemma 4.2 For two streamlines (x,<Ti(x)) and (x,<J2(x)) with 
(Ji(x) < <7 2 (x), let 

E = (0,?r) x ((7i(x),(T2(a:)). 
According to the Divergence theorem we have 

- / 1 dA = / V- V^dA = [ V^-^^ds- [ Vtp-^-^ds 

7s 7e y«n |W>| ^ 2 |W| 

= jT (jWOWx))! v/l + ai 2 (x) - |V^(x,(T 2 (x))| V / l + ^ 2 (x)) rfx 
(|^(x,<7i(a;))| (l + a' 1 2 (x)) - \^ y (x,a 2 (x))\ (l + a' 2 2 (x))) dx. 



o 



This implies that for any nontrivial streamline (x,a(x)), 

PTT PIT P (y (^£C~^ 

J \ip y (x, cr(x))\ (l + a' 2 (x)^j dx = ctt + J J ^dA, 

in view of that J n ip y (x, —d) dx = —en by the normalization (2.5). The lemma 
follows. 



Proof of Theorem 4.1 

i) In the physical variables, (X(t), Y(t)), a closed physical trajectory implies 

(X(nT),Y(nT)) = (X(0),Y(0)), (X(nT), Y(nT)) = (X(0),Y(0)) 

for some T > and all n G Z + . In the steady variables, this means 

x(0) - x(nT) = ncT, y(nT) = y{0). 

We recall that (X, Y) = (u,v), so that by periodicity - and since & = 
and v = only for x = nn - this implies 

cT = 2vrn, for some n G Z + . (4.2) 

By periodicity r = T/n is the time it takes a trajectorty x(t) in the 
steady variables to pass from x = ti to x = —ir. From (4.2) we infer 
that a physical particle trajectory is closed if and only if r = 2n/c. We 
also see from this reasoning that if r > 2n/c, then the physical difference 
X(T) — X(0) > when Y(T) = Y(0), so that the particle displays a mean 
forward drift, and contrariwise. In our case 



rir elf 

r/2 = t(0) -t(n) = - —dx 

Jo dx 

r n dx r n dx r n dx 

Jo x(x,a(x)) Jo c — u(x,a(x)) Jo \tp y (x,a(x))\ 
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According to Holders inequality 

7r 2 = ( [ dx] < [ \il> y (x,a(x))\dx [ , dX , ,, < en [ dX 

\Jo J Jo yx Jo \ipy(x,cr(x))\ Jo |^„(a:,<T(a:))| 

so that 

T > . 

c 

) For two streamlines (x,ai(x)) and (x,a 2 (x)) with a±(x) < a 2 (x) we are 
interested in the difference 



JO 



]4> y (x,(j2(x))\ \%jj y (x,ai(x))\ : 



dx. 



the positivity of which we want to prove. Let p\ and p 2 correspond to o\ 
and 02, respectively. Then we need to prove exactly that 

/ (h p (q,p 2 )-h p (q,p 1 ))dq= f h pp dqdp>Q, 
JO J s 

for E = [0,7r] x [pi,p2]- At any point either h pp < 0, so that 

fopp hqhpp ^p^qq 2h q hph q p ^ ^p^qq ^hqhphqp, 

or /ip P > 0, in which case 

Here we have used the fact that \h q \ < 1 by assumption. It thus suffices to 
investigate the sign of J- E (—hph qq + 2h q h p h qp ) dqdp. By partial integration 
in the g-variable 

' h 2 
s 



h p h qq dqdp = 2 / h q h p h qp dq dp, 



since h q = on the vertical sides of E. Hence, relying upon Lemma 3.1 i) 
and 3.2 ii), we have 

J h pp dqdp > 2 J h q h p h qp dq dp > 0. 

) Consider the quotient 

{I \1>V(X,CT 2 (X))\ ' J l ' (<7j(l) " ffl(l)) ' 

the sign of which determines the change of r and thus of the mean drift. 
The Lebesgue dominated convergence theorem can be applied to consider 
the limit a 2 — > o\, being 

d_ dx = r* ^ yy {x,a(x)) dx 

da Jo \ip y (x,a(x))\ Jo ifj^(x,a(x)) 



12 



Since ip vy = —7 > at the bifurcation point, it follows by continuity that 
the expression in (4.3) is positive in a neighbourhood of the trivial flow 
from which the nontrivial waves bifurcate. 



5 The particle trajectories 

We are now ready to discuss what our results mean for the streamlines and 
particle trajectories. We shall do so with the aid of two examples. 

5.1 Irrotational waves 

Irrotational waves display extraordinary regular features (see Figure 1). From 
bottom and up the angles between the streamlines and the horizontal plane is 
pointwise increasing, and for small enough waves this is true also for the verti- 
cal velocity. So is the maximal horizontal velocity, which for every streamline 
is attained below the crest, wherefrom it strictly decreases towards the trough. 
The surface is bounded below by a concave parabola, the curvature of which 
is determined only by the strength of gravity and the maximal horizontal 
velocity. 

In the language of particle paths, every particle traverses a non-closed cir- 
cular/elliptic orbit as the wave passes above. This forward drift is strictly 
increasing from bottom to surface, and so is the vertical movement of the 
particles (this follows from the steepness of the streamlines). At the flat bed 
there is no forward drift; the particles move equally much back and forth. As 
we move up trough the fluid, each particle still moves both backward (below 
the trough) and forward (below the crest). As the wave propagates above, the 
particle moves upwards starting from the time a trough passes until the next 
crest passes (see Figure 2). At the top of its orbit, just as the crest passes, the 
particle attains its maximal horizontal velocity. The movement then continues 
in a symmetric way, and the particle begins its descent with the horizontal 
speed strictly decreasing until it reaches its minimal value as the next trough 
passes. 

5.2 Waves of constant negative vorticity 

Some of the above features persist for waves of constant negative vorticity. 
The maximal steepness of the streamlines is strictly increasing from bed to 
surface, but we lack a proof of this property holding along any vertical line. At 
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Fig. 1. The streamlines and the velocity field for irrotational gravity waves of small 
amplitude. The steepness of the streamlines and the vertical velocity is pointwise 
increasing from bottom to surface. The horizontal velocity is everywhere decreasing 
along the streamlines, and it is increasing from bottom and up beneath the crest, 
whilst decreasing beneath the trough. 

the surface the horizontal velocity is now non-increasing from crest to trough, 
and so is it at the bottom. In between we do not know. If the vorticity is 
positive and the wave is small, then the horizontal velocity is decreasing along 
any horizontal line from crest to trough, while in the half-period to the right of 
the crest, the vertical velocity is pointwise increasing from bottom to surface. 

The surface is bounded below by the same parabola as are irrotational waves. 
The fluid particles still show a mean forward drift (except at the flat bed). 
The forward drift is however always strictly increasing from bed to surface in 
some neighbourhood of the bifurcation point (cf. [7]). We have not been able 
to verify the circular/elliptic shape of the trajectories troughout the fluid since 
there is no control of the horizontal velocity in the interior of the fluid (expect 
below the crest). At present the possibility of particles moving constantly 
forward cannot be ruled out. 



6 Some comments 

Many of the results here obtained can be extended to deep-water waves and 
solitary waves. However, it should be noted that there are important differ- 
ences between these types of waves. Notably, the investigation [5] shows that 
the particle trajectories within irrotational solitary waves differ in fundamen- 
tal ways from those in periodic waves, and in [11] it is proved that the class of 
vorticities allowed for in deep-water waves is much more restrictive than for 
finite depth. 

The detailed information presented concerning irrotational waves is in almost 
all verifiable ways consistent with the behaviour of linear waves. This is ex- 
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Fig. 2. Left: For irrotational waves and small waves of negative vorticity the forward 
drift is strictly increasing from bottom to surface. Center: To an observer stand- 
ing still as the waves passes the particles in irrotational waves traverse non-closed 
elliptic orbits, corresponding to what an observer travelling along with the wave 
understands as a streamline. Not only the the forward drift, but also the vertical 
size of these orbits is increasing from bottom to surface, and the horizontal velocity 
of each particle is strictly increasing from the bottom of the orbit to its top. 

pected but far from obvious, since a priori there is nothing to tell what "close- 
ness" between the exact wave and its linear approximation means in terms of 
sharp properties. 

In the experimental study [22] there is a picture of the particle paths within 
a progressive water wave. This shows very clearly two of the properties here 
proved mathematically: the slight forward drift and the fact that the vertical 
size of the orbits is decreasing with depth. Since all recognizable features in 
the experimental picture are consistent with our mathematical findings, it is 
likely that some of the more detailed properties we have established are also 
valid for real water waves. Hopefully, this mathematical theory can inspire the 
examination of such properties in laboratory settings. 
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